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Epidemic spreading on undirected and directed scale-free networks with correlations 
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Many complex networks have a common structure with power-law degree distributions, however 
the details such as the degree-degree correlations are different in social, technological, and biological 
systems. We numerically investigate the epidemic spreading on the network with a variety of 
correlations: the assortative, uncorrelated, and disassortative mixings. In a simulation for the mean- 
field-like susceptible-infected-recovered model, we observe different epidemic behavior according to 
the types of correlations particularly in a directed network. Our results suggest that the assortative 
connections between nodes with similar degrees enhance the epidemic spreading more significantly 
than the uncorrelated and disassortative connections between cooperative nodes with high and low 
degrees. 
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Self-organized complex networks have attracted a great 
attention to statistical physicists, computer scientists, 
and mathematical biologists, since many empirical stud- 
ies have revealed the fact that a structure is commonly 
found in social, technological, and biological networks. 
The structure is called scale-free (SF) [2] and follows a 
power-law distribution P{k) ~ k~'', 2 < 7 < 3, for the 
number of nodes with degree k. Since the heterogeneous 
characteristics of the SF network is crucial for the robust- 
ness of connectivity against failures [1] , the efficiency of 
information delivery [3], and the spread of epidemic dis- 
ease transmitted by means of social or sexual contacts, 
e-mails, Internet, and so forth [4], we expect that the 
evolutional mechanisms and the structural properties of 
the SF network [1] are useful for improving the efficiency 
and the robustness of power supply, communication, and 
economy systems. 

On the other hand, recent studies classify networks ac- 
cording to quantities of connectivity correlations of nodes 
with their neighbors [1] ; social networks tend to have as- 
sortative connections between peers with similar degrees 
[G] [7] , while technological or biological networks tend to 
have disassortative ones between those nodes with high 
degrees, namely hubs, and those with low degrees [•^>][S]. 
The properties for the epidemic incidence [!)] and for the 
percolation [10] have been compared in the considered 
forms of correlations defined as a weighted combination 
of the uncorrelated (or arbitrary correlated) term and 
the fully assortative term: a fully assortative connection 
allows only two nodes with the same degree to be con- 
nected. However the relation of the results stated above 
to the evolutional mechanisms of network is still unclear. 
Only a few analytical forms of correlations have been 
derived from the tree model [11] and the configuration 
model restricted such that at most one link exists be- 
tween any pair of nodes [12]. Although the configuration 
model produces disassortative mixing, a non-trivial dis- 
tribution of the desired degrees must be given for wiring 
in advance. Also, difficulties in estimating the condi- 
tional probability of degree-degree connectivity even from 



empirical data has been pointed out [13]. Apart from 
the evolutional mechanisms, a numerical simulation has 
shown that the introduction of assortative hub-hub con- 
nections between different local areas on a lattice has the 
effects to shorten the average delivery time and to enlarge 
the spread of infection [14]. However, neither theoretical 
solutions nor numerical studies for the epidemic behavior 
have been reported except for the above special forms of 
correlations. 

We study epidemic spreading on the SF network with 
a variety of assortative, disassortative, and uncorrelated 
connections, which are not specialized in social, tech- 
nological, or biological systems. First, we review two 
growing network models: one is called the duplication- 
divergence model [1")], which is equipped with a control 
parameter of connectivity correlations between the as- 
sortative and the disassortative mixings, and the other 
is called the directed growing model, in which the ex- 
istence of correlations is suggested [13]. Then, we will 
estimate the conditional probability of degree-degree cor- 
relations from the average realizations of the growing net- 
work models. This estimation method can be applied to 
other network models. By using the estimated probabil- 
ity, we numerically investigate the epidemic behavior for 
the mean-field-like equations of the susceptible-infected- 
recovered/removed (SIR) model. We find different epi- 
demic behavior according to the types of correlations. 

Let us consider the growing network models shown in 
Figs. 1 (a) and (b). The following procedures are re- 
peated until the network reaches to the required size N. 

Dup: duplication-divergence model [15] 

1. At each time step, a new node i' is added to the 
network. 

2. Simultaneously, a node i is randomly chosen, and 
new (undirected) links between all the neighbors j 
of i and the new node i' are duplicated. 

3. With probability qy, a link between i and i' is es- 
tablished (self- interaction). 
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FIG. 1: Growing network models, (a) In the duplication- 
divergence model, duplication links (red solid lines) between 
a new node (red circle) and all the neighbors (green circle) 
of a randomly chosen node (blue circle) , and self-interaction 
(dashed line) are generated, but some links (magenta dashed 
line) are removed, (b) In the directed growing model, new 
links from randomly chosen sources to preferentially attached 
terminals emerge. 

4. In the divergence process, each duplicated link is 
removed with probability 1 — Qe. 

These local rules are biologically plausible [10]. Note that 
larger enhances the assortativity of network generated 
by the above rules because the self-interaction means con- 
necting a pair of nodes with similar degrees. In other 
words, qy is a control parameter of the correlation. 
Dir: directed growing model [13] 

1. At each time step, a new node is added and con- 
nected from a randomly chosen node. 

2. Simultaneously, m' new links emerge from ran- 
domly chosen nodes in the network. 

3. The terminals of the new links become attached to 
nodes chosen with shifted linear preference [11]: a 
node with in-degrec k is chosen as the terminal of 
a new link with probability proportional to k + w, 
where w is a positive constant. 

The generation of new links includes the wiring between 
old nodes at each time step. In addition, to keep the con- 
nectivity of network, the first procedure is modified from 
the probabilistic addition of new nodes in Ref. [13]. If a 
multi-link between already connected nodes or self-loop 
is created, it is skipped in the directed growing model, 
while there is no such link in the duplication-divergence 
model. In a sense of reality, the sender and the receiver 
of transmitted information or objects are distinguished 
from each other on a directed link. 
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7.33 
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142 
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2.1 
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TABLE I: A set of parameters for the duplication-divergence 
and the directed growing models (denoted by Dup and Dir). 
Unnecessary parameters are marked by hyphens. The average 
degree of all nodes, the minimum and maximum degrees (in- 
degrees for directed links) of a hub node are measured over 
the 100 realizations of iV = 10^ 



In the directed growing model, the rate equations for 
the in-dcgrce distribution arc written as 

dNi m' 

= —— [(k - 1 + w)Nk-i -{k + w)Nk] + Sk.i, 
dr m + w 

(1) 

where m = m'-fl, 7 = 1 + (m + w)/m' , Sk^i is Kro- 
necker's delta, the number of nodes with degree k is de- 
noted by iVfc(r) ~ rife X T as similar to Ref. [11], and we 
obtain the solution 

^ ik-l+w)n,^^ ^ 

k -1 + W + 1 + (m + w) jm! ^ ' 

The second column of Table I shows the values of the 
parameters used in our simulation. Ass, Unc, and Dis 
denote the assortative, uncorrelated, and disassortative 
networks, respectively. We regulate these values so that 
they produce similar average degrees < fc >, because it 
is obvious that the epidemic spreading becomes larger as 
the degrees increase. Fig. 2 shows distributions of the 

def 

degree = Nk/N and the connectivity correlation 

< knn where is the conditional 

probability of the connections between nodes with de- 
grees I and k in Dup, or the connections from nodes with 
in-degree I to those with k in Dir. The degree distribu- 
tions exhibit the power-law behavior: the exponents vary 
as corresponding to the parameters 7 of Ass and Unc in 
Dir (Fig. 2 (b)), while they are close to each other in Dup 
(Fig. 2 (a)). The correlations are controlled between Ass 
and Dis in Dup (Fig. 2 (c)), while they are controlled be- 
tween Ass and Unc that is nearly uncorrelated with very 
weak correlations in Dir (Fig. 2 (d)). If the connections 
are replaced by reciprocal ones after the configuration 
of Dir, the degree distribution becomes exponential, and 
the correlation disappears. 

Next, we consider the SIR model in which individual 
nodes have three states: susceptible, infected, and re- 
covered/removed. The densities of the nodes with de- 
gree (or in-degree in Dir) k in the respective states are 
denoted by Skit) = Skit)/Nk{t), pk{t) = Ik{t)/Nk{t), 
rk{t) = Rk{t)/Nk{t). By definition, the normalization 
condition Sk{t) + Pk{t) + f'k{t) = 1 holds at each time 
t. Note that we distinguish the time-scales t and t for 
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FIG. 2: Distributions of degree and connectivity correlation, 
(a) P{k) for tiie degree in Dup, (b) P{k) for the in-degree 
in Dir, (c) < knn > for the degree in Dup, (d) < fc„„ > for 
the in-degree in Dir. The guide lines show power-law behav- 
ior of exponents 2.1 (cyan) and 3.0 (magenta). In (b), the 
yellow and black lines for the analytical solutions of Eq (2) 
are well fitting the cases of 7 = 2.1 (green) and 3.0 (red). In 
(c) and (d), the red, green, and blue lines clearly show the 
correlations corresponding to the cases of Ass, Unc, and Dis 
in Table I. The observations for very large degrees are statis- 
tically dropping and fluctuating. These are the averages over 
the 100 realizations. 



the evolution of network and the spreading of viruses. 
After the construction of a network for the Dir or Dup 
model, from an initial infected node, the epidemics is 
propagated by contacts between infected and suscepti- 
ble individuals (from infected nodes to susceptible nodes 
through directed links in Dir) at the rate b. The infected 
node is removed at the rate 6. Once an individual gets 
infected and then recovers or removed, the state is never 
changed any more. The microscopic stochastic simula- 
tion needs very expensive computation for studying the 
epidemic properties, therefore a macroscopic mean- field 
approximation is useful. Wc consider the spreading on 
the averages of randomly generated networks for each of 
Dup and Dir. 

Following Ref. [17], the mean- field- like rate equations 
for the evolution of densities can be expressed as 



dsk{t) 



dt 



-bksk{t)ek{t), 



dt 



^-6pk{t) + bksk{t)Qk{t), 



(3) 



(4) 



drkjt) 
dt 



Spkit), 



where Qk in the r.h.s of Eqs. (3) and (4) denotes the 
expectation of infection at degree fc. 



Qkit) = 



dJlf J Eti -l^PmPiit) for Dup, 



EtiPil\k)Pi{t) for Dir, 



(5) 



the factor {I — l)/l is added taking account of the fact 
that one of the links is not available for transmitting the 
infection because it was already used [17]. As applied to 
Fig. 2 (c) or (d), the P{l\k) is estimated from the aver- 
age over 100 realizations for each parameter specification 
presented in Table I. Since some degrees are missing in 
the range Kmin < k < K^ax over the realizations, we 
apply the cut-off value kc defined by Kmin- 

We numerically investigate the dynamic behavior de- 
pends on the connectivity correlations by using the 4-th 
order Runge-Kutta method for Eqs. (3) and (4) with a 
step width At = 10"^. In our simulation, we assume 
that an initial infection source is only on a hub with the 
degree k^- Hence, for the other degrees k ^ k^. wc set 
Sfc(O) = 1 and Pfc(O) = rfc(O) = 0. This assumption is 
natural since the hub is much more vulnerable against 
infection through the active communications with the 
outside. 

Figs. 3 (a) and (b) show the fraction of the infected 
nodes I{t)/N = Ik/N, and we see in the case of Ass 
(red line) it is slightly more insistent (Fig. 3 (a) for Dup) 
or larger (Fig. 3 (b) for Dir) than that in the case of Dis 
(blue line) and Unc (green line). The densities of the in- 
fected nodes pk{t) in Figs. 3 (c) and (d) show a property 
similar to that of I{t)/N. The value of pk{t) becomes 
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(d) 

FIG. 3: Time evolutions of (a) the fraction of infected nodes 
I{t)/N in Dup, (b) I{t)/N in Dir, (c) each pk{t) for k = 
1,5,40 in Dup, and (d) pk(t) in Dir at 6 = 0.1 and 5 = 0.1. 
Inset figures show similar behavior at b = 0.3 and 5 = 0.5 as 
some other example. Note that the value of Ik{t) = Nkpkit) 
is in the reverse order to k because of oc k~~' . 



larger as the degree k increases, therefore the infection 
from the nodes with high degrees are dominant in the Qk 
of Eq. (5) in particular through assortative connections. 

Fig. 4 shows the epidemic incidence R{T) =^ J2k ^k{T), 
which denotes the number of transitions from the total 
infected nodes. We set T = 100 taking the convergent 
time in Fig. 3 into account. The incidence is higher 
in the assortative networks (red line) than that in un- 
correlated (green line) or disassortative networks (blue 
line), except at 6 = 0.1 and S = 0.7,0.9 in Dup. The 
exceptions suggest that the immune rate S is an impor- 
tant factor for small b. In other words, if we set 6' = 1 
by a variable change of time t for Eqs. (3) and (4), we 
must use b' = b/S, which may be larger than 1, instead 
of < 6 < 1. These differences according to the types 
of correlations remarkably appear in Dir with smaller b. 
Furthermore, similar results are also obtained in the di- 
rected growing model based on preferentially attached 
terminals and sources by the in- and out-degrecs, respec- 
tively. 
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FIG. 4: The epidemic incidence R{T) as a function of the 
immune rate 5 in (a) Dup, and (b) Dir. The red, green, and 
blue lines are corresponding to the cases of Ass, Unc, and Dis 
at 6 = 0.1, 0.3, 0.5, 0.7 marked by the plus, cross, asterisk, and 
open square, respectively. 
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111 summary, we have numerically found different 
spreading properties on SF networks according to the 
connectivity correlations estimated from the averages of 
the growing models [13] [15]. The differences are remark- 
able in the directed model. The results suggest that 
assortative connections enhance the epidemic spreading, 
and also they could contribute to improve the efficiency of 



information delivery. In contrast, the disconnections or 
setting gatekeepers between nodes with similar degrees 
becomes a local defense strategy such as acquaintance 
immunization [18] for the spreading of viruses. Further 
works will be carefully considered for comprehending the 
effects of correlations on the spreading. 
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